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Structure determination of the indium-induced Si(lll)-(4x 1) reconstruction by 

surface x-ray diffraction 

O. Bunk,0 G. Falkenberg, J.H. Zeysing, L. Lottermoser, and R.L. Johnson 
//. Institut fiir Experimentalphysik, Universitdt Hamburg, Lumper Chaussee 149, 
D-22761 Hamburg, Germany 

M. Nielsen, F. Berg-Rasmussen, J. Baker, and R. Feidenhans'l 
Condensed Matter Physics and Chemistry Department, 
Ris0 National Laboratory, DK-4000 Roskilde, Denmark 
(21 October 1998) 

A detailed structural model for the indium-induced Si(lll)-(4x 1) surface reconstruction has been 
determined by analyzing an extensive set of x-ray-diffraction data recorded with monochromatic {Tiu: 
= 9.1 keV) synchrotron radiation. The reconstruction is quasi-one-dimensional. The main features 
in the structure are chains of silicon atoms alternating with zigzag chains of indium atoms on top of 
an essentially unperturbed silicon lattice. The indium coverage corresponds to one monolayer. The 
structural model consistently explains all previously published experimental data. 
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Considerable interest has focused recently on adsor- 
bate-induced modification of semiconductor surfaces as 
a technique to create nanoscale quantum structures of 
high perfection, fn this paper we report the forma- 
tion of quasi-one-dimensional (Id) chains on the (4x1)- 
reconstructed Si(lll) surface induced by the adsorption 
of indium and present the atomic structure as determined 
by surface x-ray diffraction (SXRD). 

Several surface reconstructions are induced by indium 
on Si(lll) (see, e.g., Kraft et at, Ref. |l|) depending on 
the coverage; the Si(lll)-(4x 1)-In reconstruction marks 
the borderline between the semiconducting, low indium 
coverage, and metallic, high-coverage, phases.tl Despite 
the fact that the indium- induced Si(lll)-(4x 1) recon- 
struction was reoorted for the first time by Lander and 
Morrison in 19650 and has been investigated with a va- 
riety of techniques subsequently, no definitive structural 
model could be established to date and some of the ex- 
perimental results are apparently inconsistent. 

In directH and inverse! photoemission investigations 
on single-domain samples and recently in inverse pho^ 
toemission investigations on a three-domain sampleEl 
the Si(lll)-(4x 1)-In reconstruction showed a quasi-ld 
metallic behavior. Additionally to this interesting elec- 
tronic property the surface exhibits an image state that 
also showed quasi-ld behavior in its strongly anisotropic 
dispersion: Along a certain direction the dispersion is 
very well described by a free electron parabola whereas 
in the perpendicular direction, the dispersion falls be- 
low the free electron parabola.B Scanning tunneling mi- 
croscopy (STMl mvestigations of the Si(lll)-(4x 1)-In 
reconstructioru'D^llil resolved zig-zag chains running in 
the (110) directions in the filled-state images and linear 
chains in the empty-state images. Tunneling data was 
acquired at bias voltages down to 0.04 V consistent with 
metallic behavior in agreement with scanning tunneling 



spectroscopy results. Adsorbed hydrogen was found to 
displace the indium atoms.uO Filled-state STM images 
of the hydrogenated substrate show a (4x1) reconstruc- 
tion with straight chains instead of the j-broad zigzag 
chains on the indium-terminated surfaccBilil From these 
obseijyations and results found in the literature Saranin 
et aln proposed a structural model. This rrpdel as well 
as the model derived by CoUazo-Davila et alS3 by apply- 
ing direct methods to transmission electron diffraction 
(TED) data are at variance with the x-ray pliatoelectron 
spectroscopy (XPS) results of Abukawa et aZJlj Since the 
core-level spectra did not show peaks corresponding to 
surface silicon atoms a complex substrate reconstruction 
can be ruled out. 

All the inconsistencies between the earlier experimen- 
tal results are eliminated with our structural model. The 
determination of the geometric structure lays the foun- 
dation for a theoretical investigation of the interesting 
electronic structure of this system. 

We used the well established experimental technique 
SXRD to determine the structure of the Si(lll)-(4x 1)- 
In reconstruction. To minimize the uncertainties induced 
by not well controlled sample preparation conditions we 
employed the unique combination of an ultra high vac- 
uum (UHV) system equipped with standard techniques 
like reflection high energy electron diffraction and low 
energy electron diffraction (LEED) facilities close to the 
wiggler beamline BW2 at HASYLAB (Hamburger Syn- 
chrotronstrahlungslabor) and a portable UHV chamber 
for the SXRD measurements. 

We used p-type (B doped) Si(lll) wafers (Wacker 
Chemie) with a resistivity of 7 f2cm and nominally no 
miscut. STM measurements on the clean Si(lll)-(7x7) 
surface showed domain sizes of more than 10000 A so 
the miscut was less than 0.02°. The carefully outgased 
sample was cleaned repeatedly by 'flashing' to ~1150°C 
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for 15-20 s and slow cooling from 900°C to room tem- 
perature. Indium was deposited from a Knudsen cell at 
a rate of ~0.4 ML/min on the sample at '--^500°C until 
the (/3 X\/3)R30° reconstruction was observed; at this 
stage the stacking fault of the Si(lll)-(7x7) reconstruc- 
tion was removed. Further deposition at a lower tem- 
perature of ~430°C (with a lower probability for indium 
desorption) yielded first the (/sT and finally the 

(4x1) reconstruction. The sample was transferred in a 
portable UHV chamber to the wiggler beamline BW2 at 
HASYLAB for the x-ray diffraction measurements. The 
incident monochromized x-rays with an energy of 9.1 keV 
impinged on the sample at a grazing angle of 0.5°. An 
in-plane data-set of 296 reflections with / ~ 0.07 was 
recorded by rotating the crystal about the surface nor- 
mal. The background subtracted integrated intensities 
were corrected for the Lorentz factor, polarization fac- 
tor, active sample area and the rod interception appro- 
priate for the z-axis geometry.Efl In SXRD measurements 
the fractional-order reflections which belong to one do- 
main do not overlap with the reflections which belong to 
the other two rotational domains. By comparing equiv- 
alent fractional-order reflections the areas of the three 
rotational domains were found to correspond to 25%, 
37% and 38% of the total surface area. The equiva- 
lent rescaled fractional-order reflections were averaged 
under the assumption of mirror- lines running along (112) 
and a systematic error of e = 13.7% was determined. 
The rods were scaled according to the corresponding in- 
plane intensity and therefore only one overall scale factor 
was necessary for the data analysis. In total the data- 
set consists of 550 symmetry inequivalent reflections, 61 
along two crystal truncation rods (ctrs), 337 along 11 
fractional-order rods and 152 in-plane reflections. 

In the following we use LEED coordinates with a = 

V2[101]eubic' b = l/mO]cnUc ^nd c = l/3[lll],^bic- 
The cubic coordinates are in units of the silicon lattice 
constant (5.43 A) and therefore \a\ = \b\ = 3.84 A and 
|c| = 3.14 A. The absolute values of the reciprocal coordi- 
nates including a factor of 27r are \a *\ = \b*\ = 1.89 A-i 
and \c* \ = 2.00 A~^. 

Usually the flrst step in the analysis of surface diffrac- 
tion data is to plot the Patterson function, i.e. the pair 
correlation function of the electron density. The Patter- 
son function projected in the surface plane can be cal- 
culated from the in-plane reflections which have a small 
momentum transfer perpendicular to the surface. A con- 
tour plot of this function is shown in Fig. ^ and is quali- 
tatively in good agreement with that previously reported 
by Finney et ala Each peak in the Patterson function 
corresponds to an important interatomic vector in the 
surface reconstruction. Since indium (Z^ = 49^) is a 
much stronger scatterer than silicon (Z^ = 14^) the peaks 
correspond to interatomic vectors between indium atoms. 
The three interatomic distance vectors which can be seen 
in the Patterson function in Fig. |^ indicate that at least 
three difi^erent indium sites must be involved in the recon- 



struction. It is highly unlikely that the indium coverage 
can be 0.5 ML (monolayer) since that would only cor- 
respond to two indium atoms per unit cell. It has been 
argued that a reconstruction of the substEate might lead 
to fuilher peaks in the Patterson functionta but the XPS 
datalij make this unlikely. 

During the course of the data analysis we tested sev- 
eral possible initial indium atom conflgurations and de- 
termined the substrate structure in the subsequent refine- 
ment. The model found is shown in Fig. |^. The correct- 
ness of the model is proven both by the good agreement 
of the measured data with the intensities calculated from 
the model structure shown in Fig. |^ and by the overall 
reduced value of 1.5. We will now describe the build- 
ing blocks of the model and resolve some of the apparent 
inconsistencies between the previously reported experi- 
mental results. 

The model shown in Fig. ^ consists of a zigzag chain 
of silicon atoms as found in the 7r-bonded chain raodel 
for the (2x1) reconstruction of the Si(lll) surfacetj on 
top of an essentially unreconstructed substrate. This 
model is consistent with the XPS investigations which 
showed that no strong siliccp. surface component was 
present in the Si 2p spectra.E^ Even in the absence of 
indium atoms the silicon chains on the substrate pos- 
sess a (4x1) periodicity and it is highly Ijtoljy that the 
chains observed after hydrogen adsorptionU'EJ are made 
up of silicon atoms. The Si(lll)-(4x 1)-In reconstruc- 
tion is completed by adding two zigzag rows of indium 
atoms in the space between the silicon chains. The ar- 
rangement of the indium atoms is similar to the arrange- 
ment proposed on the basis of ^-probe Auger electron 
diffraction inii£Btigations,tII and also used in a previous 
SXRD studyllj There are two inequivalent types-of in- 
dium atoms in agreement with the XPS results. E3 The 
indium atoms next to the silicon chains are probably 
covalently bonded to the silicon chain atoms. For the 
inner indium atoms the bonding configuration is not so 
obvious. In Fig. ^ a bonding configuration of covalent 
Px,y,z bonds to the neighboring indium atoms and down 
to the silicon substrate under an angle of approximately 
90° are shown. The nearest neighbor distance of the in- 
dium atoms is within the range of 2.98 A to 3.14 A. Elec- 
tron counting can not be strictly applied to determine the 
bonding configuration because the surface has Id metallic 
character. If we assume that the bonding configuration 
shown in Fig. ^jis a first approximation to the more com- 
plex real configuration then it is evident that the surface 
free energy of this reconstruction is lower than for other 
models since there are no silicon dangling bonds and all 
the group III indium atoms are trivalently bonded. The 
present SXRD data do not permit an accurate determi- 
nation of the bond charge densities; the detailed bonding 
configuration will have to be the subject of a future the- 
oretical investigation. 

An important question is how the model shown in 
Fig. ^with four indium atoms per unit-eell (1 ML) can be 
used to explain the results of the TEDI13 and impact coUi- 
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sion ion scattering spectrometry (ICISS) investigationsH 
which resuhed in models with only two indium atoms 
per unit cell. In the TED study direct methods were 
used to analyze the data and electron density maps with 
indium atom configurations compatible with the TED 
data were presented including plausible configurations 
with four indium atoms (e.g. Fig. 2 (f) in Ref. [ij). The 
1/2 ML model which was compatible with the TED 
data is at variance with the SXRD data as shown in 
Fig. H (c) where the dashed hues calculated using the 
1/2 ML model do not reproduce the measured data ad- 
equately. A preliminary data analysis of the SXRD data 
using direct methods resulted in the indium atom con- 
figuration shown ip-|Fig. ^thereby indicating the correct- 
ness of our model.Ea The TED data were not sufficient to 
include dynamical diffraction effects and the resolution 
was probably not good enough to rule out the wrong 
indium atom configuration.tS The reduced value for 
the model proposed by Saranin et ala is 7.3 and there- 
fore this model can also be definitely eliminated. The 
side view of the Si(lll)-(4x 1)-In reconstruction shown 
in Fig. |^(b) and the atomic coordinates given in Table || 
show that the atoms in the topmost layer have three dif- 
ferent heights. The highest atoms are the indium atoms 
bonded to the silicon chain, the inner indium atoms of 
the indium "stripe" are lower and the silicon chain atoms 
arc lowest. The existence of different heights is in agree- 
ment with the, results of a previously performed STM 
investigation.Eil Low energy ion scattering is very sensi- 
tive to the topmost surface layer. The configuration of 
the topmost indium atoms bonded to the silicon chains 
strongly resembles the 1/2 ML model with indium atoms 
on H3 and T4 sites that has been found to reproduce 
the ICISS data collected using 2 keV Li+ ions.S Another 
1/2 ML model with a zigzag chain of indium atoms on 
T4 sites also gave reasonable agreement with the ICISS 
dataQ since such chains are present in our model as shown 
in Fig. |. 

In summary, the structure of the Si(lll)-(4x 1)-In re- 
construction has been determined using surface x-ray 
diffraction. The quasi-one-dimensional character of this 
surface reconstruction is given by zigzag chains of sili- 
con atoms on top of an unreconstructed silicon substrate 
and four indium atoms per unit cell (1 ML) arranged in 
two zigzag chains in the gap between the silicon chains. 
The indium atom arrangement may also be regarded as 
being quasi hexagonal. We have shown that most of 
the previously published experimental data are consis- 
tent with the new structural model. We hope that the 
atomic coordinates given in Table || will pave the way 
for detailed theoretical investigations of the interesting 
electronic structure of this system. 
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position [LEED-coord.] 


deviation d [A] 


|d| [A] 


In 


(0.11,0.06, 0.86) 






In 


(0.86,0.93, 0.85) 






In 


(1.53,0.77, 0.99) 






In 


(3.43,0.22, 0.99) 






Si 


(2.28,0.14, 0.73) 






Si 


(2.71,0.86, 0.76) 






Si 


(0 31 65 -0 26"! 


(M 1 1 -D 05 -0 03") 

1 l^.-L-Lj \J.\J\J-^ \J-\JtJI 


0.10 


Si 


^3 96 98 -0 00"! 


(M 15 _n 08 -0 01 "> 


0.13 


Si 


(1.30,0.65,-0.25) 


(-0.12,-0.06, 0.01) 


0.10 


Si 


(0.96,0.98, 0.04) 


(-0.16,-0.08, 0.13) 


0.19 


Si 


(2.29,0.64,-0.33) 


(-0.17,-0.09,-0.24) 


0.28 


Si 


(2.01,0.00,-0.03) 


( 0.03, 0.02,-0.10) 


0.10 


Si 


(3.29,0.65,-0.23) 


(-0.16,-0.08, 0.05) 


0.15 


Si 


(2.95,0.97,-0.01) 


(-0.20,-0.10,-0.03) 


0.18 


Si 


fO 66 33 -1 24"! 


(-0 04 -0 02 04 "1 

\ KJ.yj^^ L/.WiJ^ yJ.yj^J 


0.07 


Si 


(0 33 67 -1 01 1 
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1 yj-yjA.^ yj.yjA.^ i 


0.02 


Si 


(1.65,0.32,-1.26) 


(-0.08,-0.04,-0.03) 


0.07 


Si 


(1.32,0.66,-0.99) 


(-0.06,-0.03, 0.03) 


0.06 


Si 


(2.67,0.33,-1.27) 


( 0.00, 0.00,-0.05) 


0.05 


Si 


(2.32,0.66,-1.04) 


(-0.05,-0.03,-0.13) 


0.14 


Si 


(3.66,0.33,-1.25) 


(-0.04,-0.02, 0.00) 


0.03 


Si 


(3.31,0.65,-0.98) 


(-0.09,-0.05, 0.04) 


0.09 



TABLE I. The atom positions in the Si(lll)-(4x 1)-In re- 
construction derived from the analysis of the SXRD data. 
The positions of the atoms are given in LEED coordinates, 
the deviations from the bulk-like positions and the absolute 
values of these deviations are given in A. For the sihcon atoms 
an isotropic Debye-Waller factor of 0.5 A^ and for the indium 
atoms of 4.5 A^ was used in the data analysis. 



FIG. 1. Patterson function of the electron density pro- 
jected in the surface plane calculated from the fractional-order 
in-plane reflections. The axes are scaled in LEED coordinates 
(1.0 corresponds to 3.84 A). The distance vectors shown are 
1. (0.69,0.85), 2. (1.42,0.71) and 3. (2.00,0.50). The dashed 
line from (0,0) to (2,1) indicates the mirror line (along (112) 
in bulk coordinates). 

FIG. 2. BaU and stick model of the Si(lll)-(4x 1)-In re- 
construction in top (a) and side (b) views. Indium atoms are 
drawn dark grey, silicon atoms are drawn light grey. The con- 
tributions to the peaks of the Patterson function are shown 
by arrows. The standard LEED (4x1) unit cell is indicated 
by a dashed line. The dashed-dotted line along [112] indicates 
a mirror line. 

FIG. 3. Measured and calculated SXRD intensities: (a) 
in-plane with I = 0.07. The radii of the filled (empty) 
semi-circles are proportional to the measured (calculated) in- 
tensities. Hatched circles are scaled with a factor of 0.5. (b) 
Fractional-order rod-scans. The solid line is calculated using 
the model shown in Fig. ^with the coordinates given in Table 
|. (c) Integer order rods. The dashed lines correspond to the 
intensities calculated using the best fit to the model proposed 
in Ref. |l|. 
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